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ABSTRACT 

We conduct a deep mid-infrared census of nine massive galaxy clusters at (0 < z < 1.3) with a total 
of ~ 1500 spectroscopically confirmed member galaxies using Spitzer /IJiAG photometry and established 
mid-infrared color selection techniques. Of the 949 cluster galaxies that are detected in at least three of 
the four IRAC channels at the > 3cr level, we identify 12 that host mid-infrared selected active galactic 
nuclei (IR-AGN). To compare the IR-AGN across our redshift range, we define two complete samples of 
cluster galaxies: (1) optically-selected members with rest-frame Vab magnitude < —21.5 and (2) mid-IR 
selected members brighter than (Mg g-)-0.5), i.e. essentially a stellar mass cut. In both samples, we 
measure /ir-agn ^ 1% with a strong upper limit of ~ 3% at z < 1. This uniformly low IR-AGN fraction 
at z < 1 is surprising given the fraction of 24/im sources in the same galaxy clusters is observed to increase 
by about a factor of four from z ^ to z ^ 1 ; this indicates that most of the detected 24/im flux is due to 
star formation. Only in our single galaxy cluster at z = 1.24 is the IR-AGN fraction measurably higher 
at ~ 15% (all members; ~ 70% for late- types only). In agreement with recent studies, we find the cluster 
IR-AGN are predominantly hosted by late-type galaxies with blue optical colors, i.e. members with 
recent/ongoing star formation. The four brightest IR-AGN are also X-ray sources; these IR-I- X-ray AGN 
all lie outside the cluster core {Rproj ^ 0.5 Mpc) and are hosted by highly morphologically disturbed 
members. Although our sample is limited, our results suggest that /ir-agn in massive galaxy clusters 
is not strongly correlated with star formation at z < 1, and that IR-AGN have a more prominent role at 

Subject headings: galaxies: active - galaxies: clusters: general - infrared: galaxies - galaxies: evolution 



1. INTRODUCTION 

Our understanding of galaxy evolution in dense clus- 
ter environments continues to improve as observations 
broaden in wavelength and redshift range. One funda- 
mental observation is that dense environments at low 
redshifts are dominated by passive, early-type galax- 
ies that define a narrow red sequ ence in an optical 
color-magnitude d i agram (ISan dagc fc VisvanathanI 119781 : 
IBower et al] [T99l IHogg et all 12004), and studies of the 
cluster color-magnitude relation (CMR) show that the lu- 
minous red sequence members in galaxy cl usters have no t 
evolved significantly since z ~ 0.8 ( Rudnick et al.|[2009D . 
In contrast, the less massive cluster members continue to 
migrate to the red sequence as studies at z > 0.6 show 
there are fewer faint red galaxies in cl usters compared the 
field and to lowe r redshift clu sters (De L ucia et al.l 120071 : 
IStott et all [2007t [Rud^lick et a l. 2009). Recent observa- 
tions at z > 1.4 hav e now even found massive, st ar forming 
galaxies in clusters (ITran et al.| | 2010l: iHlIton et al. 201 0,). 

The question then remains as to what halts star forma- 
tion in cluster galaxies? Possib le enviro nmental processes 
include ram-pressure stripping (jGunn fc Gott 1972), tidal 
effect s from the cluster potential (Farouki & Shapiro 
I198ir ). and galaxy-galaxy interactions (Richstonc 1976), 
but none are completely effective at reproducing the star 
formation histories and scaling relations observed in galaxy 
clusters. One model that has proven successful is in- 
cluding feedback from active galactic nuclei (AGN) . Semi- 
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analytic models that include AGN ar e able to reproduce 
observed mass /luminosity functions (|Croton et al.l [20061 : 
Bower et al. 2006; Lagos ct al. 2008). AGN can also af- 
fect the intra-cluster medium (ICM) where models find 
that including AGN produces much better ag reement with 
observations of X-ray properties of the ICM (IBower et al.l 
[200llPuchwein et al.l[200a IMcCarthv"eFan [20To[r 

AGN feedback seems to be the ideal solution to re- 
solve many outstanding discrepan cies between models 
and observatio ns (Gabor et al. 2 0T0l : iFbntanot et al.ll2010l : 
ITevssier et al.ll201Q) : however, there is not yet clear ob- 
servational evidence that AGN contribute significantly to 
the quenching of star formation in cluster galaxies, mak- 
ing clusters into resting homes of passive galaxies. Sev- 
eral groups using primarily X-ray observations find the 
fractional density of X-ray selected active galactic nu- 
clei (X-r ay AGN) in cluster environments increas es with 
redshift ([Eastman et all [20071 : [Martini et al.l 120091 ). Simi- 
larly, using X-ray, infr ared and radio selection at z < 1.5 
IGalametz et al.l (|2009l ) find that the AGN surface-density 
(A^/arcmin^) is greater for clusters than in the field and 
that the AGN volume-density (N/Mpc^) for clusters in- 
creases with redshift. In contrast, a study of CL 0023-1-04, 
a large scale system of four galaxy groups merging at 
z ^ 0.83, do es not find an excess of X-ray sources relative 
to the field ([Kocevski et al.i r 2009l ). In general, studies are 
hampered by the small number of X-ray AGN and the need 
to isolate a large sample of cluster galaxies at higher red- 
shifts. Another important issue is that differ ent diagnos- 
tics select different populations of AGN (.Hart et al.l 120091 : 
IHickox et an[2009l : [Griffith fc Stern|[20T0[ ) and so no single 
approach will be complete. Thus whil e CL 0023 -I- 04 does 
not have an excess of X-ray sources, ILubin et al.l ([20091 ) 
does find a population of passive (no detectable Ha) mem- 
bers with [Oil] emission that may be due to AGN. 
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TABLE 1 
Cluster Properties 



Cluster 


RA Dec'' 


2-Range'' 




A^IRAC'^ 


logiLxT 


Status* 


Coma 


12 59 35.7 


+27 57 34 


(0.013 - 


0.033) 


348 


262 


43.0 


Relaxed 


Abell 1689 


13 11 29.5 


-01 20 17 


(0.17 


- 0.22) 


81 


73 


43.3 


Relaxed 


MS 1358+62 


13 59 50.4 


+62 31 03 


(0.315 - 


0.342) 


171 


133 


43.0 


Unrelaxed 


CL 0024+17 


00 26 35.7 


+17 09 43 


(0.373 - 


0.402) 


205 


75 


42.5 


Unrelaxed 


MS 0451-03 


04 54 10.9 


-03 01 07 


(0.52 


- 0.56) 


242 


90 


43.3 


Relaxed 


MS 2053-04 


20 56 21.3 


-04 37 51 


(0.57 


- 0.60) 


132 


87 


42.8 


Unrelaxed 


MS 1054-03 


10 57 00.0 


-03 37 36 


(0.80 


- 0.86) 


153 


120 


43.2 


Unrelaxed 


RX J0152-13 


01 52 43.9 


-13 57 19 


(0.81 


- 0.87) 


147 


80 


43.3 


Unrelaxed 


RDCS J1252-29 


12 52 54.4 


-29 27 18 


(1.22 


- 1.25) 


38 


29 


42.8 


Unrelaxed 



Note. — All clusters in this sample have Mvir > 5 X 10** M©. 
" Right ascension and Declination are for J2000. 

^ Range of re dshifts for spectroscopically confirmed members from ISaintonge et al] II2008I ) and 
IDemarco et~al] f2007i . RDCS J1252-29) 

Total number of spectroscopically confirmed cluster members. Redshifts are from Rines et al. (2003, 
Coma), Due et al. (2002, Abell 1689), Fisher et al. (1998, MS 1358), Moran et al. (2005, CL 0024), 
Moran et al. (2007a, MS 0451), Tran et al. (2005a, MS 2053), Tran et al. (2007, MS 1054), Demarco 
et al. (2005, RX J0152) and Demarco et al. (2007, RDCS J1252) respectively. 
"* Number of cluster galaxies with detections in at least three IRAC channels. 

" Bolometric ICM X-ray luminosities (ergs s"*) from Holden et al. (2007, Coma, MS 1358, MS 2053, 
MS 1054, RX J0152), Bardeau et al. (2007, A1689), Donahue et al. (1999, MS 0451), Zhang et al. 
(2005, CL 0024) and Rosati et al. (2004, RDCS J1252). 

* Dynamical state of each cluster determined from redshift distributions and X-ray & weak lensing 
profiles. Unrelaxed systems are those that show signs of a cluster scale merger. 



With the launch of the Spitzer Space Telescope, mid- 
infrared selection has become an efficient alternative 
method for identifying AGN that complements existing X- 
ray, radio, and optical search techniques. Note that mid-IR 
color selection is parti cularly effective at selectin g high- 
Eddi ngton AGN (e.g.. IHickox et al l 12009!: Hopkin s et al.l 
I2009[ ). A single stellar population (SSP) has a declining 
Rayleigh- Jeans tail that longward of f /i m makes thes e type 
of SEDs blue in the mid-infrared (Zhan g et al.ll200^ while 
an S EP with an activ e galactic nucleus has a rising power 
law (jElvis et al.lll994D and so will be red. Advantages of 
selecting with mid-infrared wavelengths include decreased 
attenuation from dust and sensitivity to AGN spanning 
a broad range in redshift (0 < z < 2) (jStern et all 120051 : 
ILacv et ani2004l ) . In the mid-IR, red colors are due to ther- 
mal emission from heated dust, either by radiation from 
stars or an active nucleus, and so dust heated by young 
stellar populations, e.g. ultra-luminous infrared galax- 
ies (UL IRGs), can contarn inate a mid-IR selected AGN 
sample (jPonlev et al.|[2008[ ). However, (U)LIRGs are ex- 
tremely rare in massive galaxy clusters with only one con- 
firmed U LIRG out of thousands of known cluster g alaxies 
at z < 1 (ISaintonge et al.ll200l IChung et aLl[20ToD . 

We expect that most cluster galaxies will populate a 
tight distribution in mid-IR color space because mas- 
sive galaxies in clusters formed the bulk of th eir stars at 
z > 2 and have evolved passively sinc e (Bower et al.1119921: 
van Dokkum et al.l [1991 ITran et al.l [2007- .Muzzin et al.l 



2008( 1. However, if AGN are essential for quenching star 



formation in cluster galaxies, the AGN fraction should 
increase with redshift because the IR luminosity func- 
tion and fraction of star forming galaxies in massive clus- 
ters increases b y a facto r > 4 from z ~ to 2; ~ 1 
(|Saintorige et al' 2008; B ai et al.l[20T0[ ). We can identify 
AGN using mid-IR colors to determine if the IR-AGN frac- 
tion evolves with redshift and, combined with known X-ray 
AGN, also study the properties of cluster AGN to obtain 
a more complete picture of AGN in galaxy clusters. Note 
that selecting AGN with optical diagnostics becomes quite 
challenging at z > 0.4 because key spectral features shift 



into the near-infrared, and extremely deep radio obser- 
vations are significantly more time intensive than mid-IR 
imaging. 

Motivated by these issues, we present the first exten- 
sive mid-infrared census of massive galaxy clusters at 
< z < 1.3 (Table [T]). Our clusters are the most massive 
systems known (Mvir ^ 5 x lO^'^M©) having comprehen- 
sive spectroscopic information for member galaxies as well 
as deep optical photometry (see §2.21) and so provide the 
most representative sample of galaxy clusters in the range 
< z < 1.3. This paper is organized as follows: in §2 
we discuss the reductions of the Spitzer mid-infrared and 
optical data for each cluster. In §3 wc describe how we 
select AGN and results for each cluster. In §4 we measure 
the IR-AGN fraction as a function of redshift and discuss 
properties of their host galaxies, and we present our con- 
clusions in §5. All magnitudes are in the AB system unless 
otherwise noted. We assume a ACDM cosmology through- 
out with 0.27 , 17a= 0.73 and h = 0.71 

2. DATA AND REDUCTIONS 

2.1. Spitzer IRAC 

Mid- infrared observations were taken with the IRAC in- 
strument (Fazi o et al.ll200^ o n board the Spitzer Space 
Telescope (Werner et al.ll20'M ) and are publicly available 
on the Spitzer archive. IRAC observes in four mid-infrared 
channels centered at 3.6, 4.5, 5.8 and 8.0 /im. These chan- 
nels have transmission functions such that all emission 
between 3.1 and 9.4 fim will be detected. Furthermore, 
the edges of the channels are steep, so that as an emis- 
sion feature becomes redshifted out of one bandpass it will 
promptly shift into the adjacent one. This ensures that 
spectral features (such as PAH emission from dusty star 
formation) can be traced fairly easily over a wide range in 
redshift. 

For most of the cluster fields (Table [1]), multiple IRAC 
observations were conducted at different times. In the in- 
terest of increasing depth we make use of all available data. 
This also frequently served to increase the area of coverage 
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allowing more cluster galaxies to be detected. 

2.1.1. Mosaicking 

Each observational campaign is composed of a series 
of dithered images referred to as Basic Calibrated Data 
(BCD). The depth of a final mosaic depends on the com- 
bination of exposure time and number of frames per sky 
position. BCDs are single-frame images that have been 
reduced and flux-calibrated on a basic level by the Spitzer 
pipeline. We performed fur ther processing and m osaicking 
using MOPEX(18.3 rev 1, IMakovoz et al|[200l . a set of 
reduction and analysis tools designed by the Spitzer Sci- 
ence Center. The software includes modules that flag and 
remove spurious detections that are not accounted for in 
the automatic pipeline reduction. Prior to mosaicking, im- 
ages were inspected for artifacts such as muxbleed and col- 
umn pull- up/down and were mitigated using the cosmetic 
module from MOPEX. 

Overlapping is another necessary step in the image re- 
duction process available as a package in MOPEX. Within 
a set of dithered frames, the individual backgrounds may 
vary to the point that the background of a final mosaic 
shows a checkered pattern/gradient. Performing over- 
lap matches the backgrounds of all input frames and so 
removes this effect. After overlapping, BCDs are then 
stacked. Weight maps were also obtained for each mosaic 
as coverage varied between campaigns. 

2.1.2. Photometry and Completeness 

Fixed-aperture photometry was carried out on each mo- 
saic using SExtractor 2.5.0 (Bertin & Arnouts 1996). In 
order to analyze cluster galaxies at varying z equally, aper- 
tures were chosen at constant proper si zes ac cording to 
cluster redshifts (except for Coma, see §3.2.ip . Ideally, 
apertures should be small so as to isolate the central en- 
gine and reduce possible contamination from, e.g. ex- 
tended star forming regions in the host galaxy. However, 
the IRAC point spread functions (PSFs) are ^^1.66, 1.72, 
1.88 and 1.98" in diameter for channels 1—4 respectively, 
and fluxes determined from apertures comparable to the 
IRAC pixel size (~ 1.22" px~^) are subject to sampling 
errors. These caveats constrain the minimum size of a 
reasonable aperture. Ultimately, we choose a radius of 
~12.6 kpc, which corresponds to an aperture diameter of 
3" for the most distant cluster (see Table [2] for all aper- 
ture sizes). We note that a 12.6 kpc aperture is much 
larger than would be "ideal" for this type of a study; how- 
ever, we feel that accepting this limitation is preferred over 
varying the physical sizes of apertures, which would likely 
introduce a selection bias. 

Aperture correc tions were determined from IRAC 
calibration s tars (jReach et al.l [2005) as discussed in 
lAshbv et al.l (|2OO90 , except for Coma where theoreti- 
cal aperture corrections for extended sources were used. 
Although some galaxies in our other low-z clusters are 
resolved, the discrepancy between extended- and point- 
source aperture corrections is not enough to affect our re- 
sults. Further discussion regarding these corrections can 
be found in the IRAC Instrument Handbook. IRAC fluxes 
are calibrated based on 24" (diameter) apertures and an 
appropriate correction needs to be applied to photometry 
using apertures of different sizes. Aperture corrections are 
defined as the difference between the magnitude of a point- 
sources from a 24" aperture and the magnitude from the 
aperture of interest. We determined aperture corrections 
from five standard stars (HD-165459, 1812095, KF06T1, 
NPMlp66.0578, NPMlp67.0536) using the average as the 



TABLE 2 
IRAC Photometry 



Cluster 


Aperture'' 


3.6/jm ^ 


4.5/xm ^ 


5.8Atm ^ 


S.OAtm b 


Coma 


15" 


20.3 


20.1 


18.8 


18.1 


1689 


7.7" 


19.6 


19.9 


19.7 


19.0 


1358 


5.3" 


20.8 


21.1 


20.4 


20.2 


0024 


4.8" 


20.6 


20.9 


20.2 


20.0 


0451 


4.0" 


20.9 


21.0 


20.3 


20.0 


2053 


3.8" 


20.8 


20.9 


20.3 


20.1 


1054 


3.3" 


21.1 


21.2 


20.7 


20.5 


0152 


3.3" 


21.1 


21.4 


20.1 


20.3 


1252 


3.0" 


21.0 


21.3 


20.6 


20.6 



^ Fixed circular aperture corresponding to physical diameter of 
~ 25 kpc at cluster redshift except for Coma where the aperture 
corresponds to 6 kpc ; see §3.2 for more details. 
^ AB magnitude corresponding to 80% completeness limit; see 
^2.1.2l for more details. 

final value. Because Galactic extinction is negligible at 
these wavelengths (< 0.01 mag) the corrections are less 
than the measurement uncertainties, we do not correct the 
IRAC magnitudes; however, we do correct for Galactic ex- 
tinction in the optical filters. 

Completeness was measured using the gallist and mkob- 
jects modules in IRAPQ. For each mosaic, 1000 artificial 
galaxies in half-magnitude bins between 16 < to < 25 were 
distributed randomly. Source extraction was then carried 
out for these fake sources with identical parameters as used 
for the real sources. Fake sources that were extracted with 
magnitudes brighter than their input value were discarded 
as blendings. In order to get a sense of the completeness 
within each cluster, sources were added in 1.5x1.5 Mpc 
boxes centered on the core of each cluster. Measurements 
of completeness are shown in Figure [T] and Table [2j 

In performing source matching, we start with the coor- 
dinates of a galaxy from optical images. We then search 
for the nearest source in each IRAC catalog within a 2" 
radius, i.e. slightly larger than the IRAC PSF. In a few 
cases, e.g. in the cores of the high redshift clusters, there 
is some blending of IRAC sources. We deblend and sepa- 
rate sources using the SExtractor parameters DEBLEND 
NTHRESH=64 (for afl channels) and DEBLEND MIN- 
CONT=0 and 0.005 (for channels 1,2 and 3,4 respectively). 
Visual inspection confirms that these parameters are ef- 
fective at deblending sources with only 2 pairs/groups of 
galaxies still blended in the 2 > 0.8 clusters; only one of 
these has an IR-AGN signature (see §3.2.8). 

2.2. Optical Photometry and Spectroscopy 

Catalogs of optical and near infrared photometry as 
well as optical spectroscopy were obtained from multiple 
sources. Observed photometry was converted to rest-frame 
values using KCorrect v O.2.1, an extension of kcorrect 
(jBlant on fc Rowcis 20071 developed for Python by Taro 
Sato. Extensive spectroscopic catalogs were used to con- 
firm membership for each cluster. General properties of 
each galaxy cluster are shown in Table 2. To summarize: 

• Coma : Photometry was taken from IMobasher et al.l 
(j2001h which consists of measurements in the B and 
R bands. Completeness was assessed to be 22.5 mag 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatory (NOAO), which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. — Completeness functions for each cluster field showing the fraction of recovered sources versus their input magnitude. For each cluster, 
1000 generated galaxies were randomly distributed in 0.5 magnitude bins in a 1.5x1.5 Mpc^ box centered near the cluster core. Sources that 
were extracted with magnitudes brighter than their input value were rejecte d as blendings. The 3.6 and 4.5/im channels exhibit a more gradual 
decline in completeness due to source confusion from greater crowding (see IBarmbv et al.ll2008t ) which increases the likelihood of blendings. 



and 21 mag in B and R respectively. Spectroscopy 
for the Coma cluster was taken as part ol the Clus- 
ter a nd Infall Nearby Survey (CAIRNS: iRines et alJ 
120031 ). Galaxies targeted in this survey were selected 
from digitized images of the POSS I 103aE (red) 
plates which are complete down to E=15-16. 



i Abell 1689: Data for Abell 1689 Puc et all [200l 
were taken as follow -up to the photom etric observa- 
tions conducted by IDve et all (|2001| ) . Photometry 
was acquired for the BVR bands and is complete to 
23.0, 22.7 and 22.7 mag respectively. Spectra were 
obtained for most(~75%) of the photometric cluster 
members at i? < 17.75 mag which drops to ^40% at 
R < 19.5 mag. 

' Mg 1358+62- Obs^ ervations of MS 1358 taken from 
IFisher et al.l () 19981 ) including photometry in the V 
and R bands and spectroscopy. Spectroscopic com- 
pleteness was determined to be >80% at i? < 21 mag 
when compared to photometric observations which 
were complete to i? ^ 23.5 mag. 



• CL 0024+17 and MS 0451-03: Data for CL 0024 
and MS 0451, including photometry from the HST 
WFPC2 i n struni ent, are discussed in d etail in 
ITreu et al.l (pOOl and IMoran et al.1 (pOOSl 2007). 
Photometry was measured to be complete to / 25 
(Vega mags)EI for CL 0024 and spectroscopic com- 
pleteness was found to be > 65% at / < 21.1 and 
/ < 22.0 mag for CL 0024 and MS 0451 respectively. 
Additional ground-based photometry was also ob- 
tained in the BVRIJKg bands reaching Sa depths 
of 27.8, 26.9, 26.6, 25.9, 21.6 & 19.7 mag for CL 0024 
and 28.1, 27.0, 27.3 25.9, & 20.2 mag for MS 0451 
(|Moran et al.l[2007bf ). 

• MS 2053-04 and MS 1054-03 : Spectroscopy for MS 
2053 is detailed in lTran et al.l (|2005a) and complete- 
ness determinations were assessed according to sam- 
pling and success rates. The success rate is defined 

2 Here the 7-band refers to the F814W filter from the WFPC2 
instrument 
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Fig. 2. — Combined IRAC color-color plots as defined by I Stern et al.l I I2005I . left) and ILacv et al.l 120041 . right) derived from 3" apertures. 
Plotted are all detected sources in the IRAC mosaics for each cluster field. Curves in both panels are redshift-tracks for M82 (a local starbursting 
galaxy; blue), VCC 1003 (a local passively evolving ga laxy; red), Mrk231 (Seyfert 1 A GN; black) and a Seyfert 2 template (pink); tracks begin 
at ^ = marked with Xs and go to ^ = 2 UDevriendt et al.lll99g| : IPolletta et al|[2007t) . The areas enclosed by the dashed lines are empirically 
defined regions designed to select galaxies dominated by emission from an AGN. Sources that fall in the wedge on the left are overlayed with 
cyan squares on the right whereas sources in the wedge on the right are overlayed with green diamonds on the left. Only ~33% of Lacy IR-AGN 
are selected as Stern IR-AGN, whereas ~89% of Stern IR-AGN are selected as Lacy IR-AGN. Furthermore, the Stern criteria determine an 
AGN fraction of 15%, whereas using the Lacy criteria it is 40 %. We adopt the Stern IR selection for our analysis as it seems to suffer less 
contamination from non-AGN sources. See lBarmbv et aT\ l|2006l ) and lDonlev et aLl H2008l 'l for a more in depth analysis of mid-IR color evolution 
in galaxies with and without AGN. 



as the number of spectroscopic redshifts obtained di- 
vided by the number of spectroscopic targets. Spec- 
troscopic completeness was determined to be ^70% 
at rn8i4 < 22 mag. Similar m ethods were applied for 
MS 1054 (ITran et al.l [20071 ) which found complete- 
ness to be >75% at nisi4 < 21.2 mag. Photometry 
for both MS 2053 and MS 1054 were acquired from 
the HST/WFPC2 F606W and F814W fihers. 

• RX J0152-13 and RD CS J 12 52-29: Photometry 
for RX J0152 (Blakcsleiliril] [2006) was obtained 
from the ACS instrument onboard the HST in the 
F625W, F775W and F850LP bandpasses. Incom- 
pleteness for these obse rvations begins to se t in 
at ~23.5, 22.5 and 22. IDemarco et all (pOOSi ) de- 
termined spectroscopic membership for RX J0152 
confirming 102 cluster galaxies out of 262 tar- 
gets. Observations of RDCS 1252 is outlined in 
IDemarco et al.l (|2007l ). Photometry was taken in the 
BV Ritj^z^^qJ K s filters reaching 5(t limiting magni- 
tudes of 26.5 & 26 mag in the J k, Ks filters re- 
spectively. The spectroscopic success rate for RDCS 
J1252 was found to show a rapid decline at Kg = 21.5 
mag, dropping from 85% to 20%. 

3. RESULTS 

3.1. IRAC Color Selection of ACN 

Mid-IR emission from an AGN is widely ac- 
cepted as thermal continuum from circumnucl ear dust 
([Andreani et al.l 120031: IKuraszkiewicz et al] 120031 ). As ra- 
diation from the accretion disk bombards the surround- 
ing dust, it is heated to temperatures in the range of 
~ 20 - 1000 K (belo w dust subhmation at T ~ 2000 
K. [Sanders et al.l [19891) . IRAC colors alone can be an ef- 
fective method for separating star forming ga laxies from 
those hosti ng AGN at redshifts up to ^ ~ 2 ([Lacv et al.l 
[2004: .Stern eTaLl [2001 . It is worth noting that > 50% 



of a galaxy's mid- IR emission must originate from the nu- 
clear component (iHlckoxet al.ll2007l: [Hopkins et al.1120091: 

'Atlee et al. '2011") in order to be selected by the criteria 
of Stern et al. (2005). Thus, due to various limitations in 
measuring IRAC fluxes in our sample, we are only able to 
select galaxies that are dominated by AGN emission in the 
mid-IR. 

The optical to mid-IR spectral energy distribution 
(SED) of AGN are typically well characterized by a rising 
power law with a few notable PAH features, thus caus- 
ing thes e objects to appear increasingly red in the IRAC 
window (IPolle tta et al.ll2007l ). Consequently, AGN tend to 
populate a separate region in mid-IR color space. Further- 
more, studies have shown that various types of Seyf erts 
have very similar mid-IR properties (Ga ndhi et al.i ,20091 
and references therein) . Radiation at these wavelengths is 
relatively insensitive to extinction and thus gives a reliable 
measure of reprocessed emission from the central engine. 
However, it is worth noting that highly obscured sources 
{Ay > 30) may be pushed outside th e IRAC selection 
wedge (see Fig. 1 of lHickox et~all [20071) . 

Figure [2] shows the em pirical co lor selection criter ia from 
both lStern et al.l (|2005l left) and ILacv etUI (pOi, right) 
for all IRAC sources detected in our nine cluster fields 
(including the cluster galaxies). Included in Figure [2| are 
color tracks of template galaxy SEDs from lDevriendt et al.l 
(jT999j): M82 (local starbursting galaxy; blue), VCC 1003 
(local passively evolving galaxy; red), Mrk231 (Seyfert 1 
AGN; black), and a Seyfert 2 template (pink). Tracks 
begin at z ~ m arked with Xs and go to z = 2 
Pollettaet al.ll2007D. W e refer the reader to lBarmbv et al.l 
20061) and lDonlev et al.l ([2008D for a more detailed discus- 
sion of the mid-IR color evolution for galaxies including 
limitations of mid-IR selection. 

T o compare the sele ction methods of ILacv et al.l ([20041 ) 
and [Stern et al.l ([20051 ) we determine IR-AGN fractions for 
all galaxies (field and cluster) detected in all four bands of 
our IRAC imaging. We obtain fractions of 40% and 15% 
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Fig. 3. — IRAC color-color plots used to select active galactic nuclei (IR-AGN: IStern et al.|[2005l) for the galaxy clusters in our study. Cluster 
redshift and number of confirmed members are shown in each panel. Data points correspond to morphologically classified early-type members 
(red circles), late-types (blue diamonds) and unclassified/mergers (black circles). Upper limits are determined for galaxies lacking detections 
in a particular channel (green and pink arrows, see [[3]|. X-ray sources are indicated as stars. Mean uncertainties are shown in the lower-right 
corner of each panel. Early- type galaxies with blue IRAC colors populate the lower-left region in each plot (the "passive cloud") whereas 
IR-AGN populate the area enclosed by the dashed lines. Our cluster IR-AGN are predominantly hosted by late-type galaxies. 



for the Lacy and Stern criteria respectively. In Figure [2l 
the green diamonds in the Stern plot correspond to sources 
that are selected as AGN using the Lacy criteria, whereas 
cyan squares in the Lacy plot correspond to sources that 
are selected as AGN using the Stern criteria. Approxi- 
mately 89% of the Stern-selected AGN are Lacy-selected 
AGN; however, the reverse shows that only ~33% of all 
Lacy-selected AGN are also selected based on the Stern 
criteria. The higher IR-AGN fraction measured using the 
Lacy criteria is not surprising given that the track for M82 
falls in the Lacy wedge. Because the Lacy criteria do not 
exclude starburst galaxies as effectively as the Stern crite- 
ria, we adopt the Stern criteria throughout the rest of this 
paper. 

3.2. Individual Clusters 

IRAC color plots for the nine massive galaxy clusters 
are shown in Figure [31 Only galaxies that have been 
spectroscopically confirmed as members with > 3cr detec- 
tions in at least 3 IRAC channels are shown. Data points 
indicate morphologically classified elliptical/SO galaxies 
(red circles), late- type galaxies (blue diamonds), unclas- 
sified/merger (black circles), sources with no detection in 
channel 4 only (pink arrows) and no detection in channel 3 
only (green arrows). For sources that lacked detections in a 



single bandpass, upper/lower limits were determined by as- 
suming the 80% completeness magnitude for the respective 
bandpass (Table [2]). Applying these limits mostly reveals 
a fainter population of passive galaxies. Considering that 
channels 1 & 2 (shorter wavelength) probe to fainter mag- 
nitudes than channels 3 & 4 (longer wavelength), galaxies 
in the "passive cloud" (with declining mid-IR SEDs) tend 
not to be detected at longer wavelengths while galaxies 
with IR-AGN (with increasing mid-IR spectra) are more 
likely to be deteced in all four channels. This is why no 
potential candidate AGN are identified by our limit deter- 
minations. We use optically-determined coordinates for 
cluster galaxies to locate their IRAC counterparts. Using 
a matching radius of 2" (6" for Coma) we find the rate of 
detecting a false positive to be < 1%. 

3.2.1. Coma 

The Coma cluster is one of the richest and most closely 
studied galaxy clusters and is known to be dominated by 
passively evolving systems with early-type morphologies 
ijMichard fc Andreonll2008l) . IRAC imaging for this cluster 
in all four channels covers roughly a 51.1' x 62.5' region cen- 
tered on NGC 4874, a field of view that includes 348 spec- 
troscopically confirmed cluster galaxies from the CAIRNS 
(iRines et al..,20G3 ). Determining photometry with Spitzer 
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data for resolved galaxies is not as straightforward as for 
the mostly unresolved galaxies in the distant clusters (see 
the IRAC Instrument Handbook for details). We decide 
against using flexible-aperture photometry from SExtrac- 
tor because of the differently sized apertures used for the 
same galaxy across the IRAC channels. Instead, we use 
a constant 15" diameter aperture (6 kpc at Coma's mean 
redshift) and apply aperture corrections to the extended 
objects as detailed in the IRAC Instrument Handbook. 

Not surprisingly, we find that an overwhelming majority 
of Coma galaxies occupy the "passive cloud" in IRAC color 
space and well below the AGN wedge. The dispersion in 
the [5.8] — [8.0] color is likely due to PAH features from star 
formation at 6.2/xm and 7.7/^m being detected in IRAC's 
8.0/im bandpass; note that Coma's proximity means we 
detect even the faintest members. There is one member 
with a significantly redder [3.6] — [4.5] color that is a disk 
galaxy viewed at an intermedi a te an gle, but it is not an 
IR-AGN from the IStern et all ()2005[ ) criteria and is not 
cla ssified as an optical A GN in the recent survey of Coma 
bv lMahaian et al.l (pOTOt ). 

3.2.2. Abell 1689 

Spectrosco pically confirmed m embers and their photom- 
etry are from lDuc efall (pOOl . The scatter in [5.8] -[8.0] 
color among members can be attributed to PAH features 
at 3.3, 6.2 and 7.7/im in dusty star forming galaxies where 
the latter two features both shift into the 8.0/im channel. 
However, the 3.3/im PAH feature shifts to the boundary 
between channels 1 & 2. We find one candidate IR-AGN 
for this low redshift cluster that has also been c lassified as 
a Sey fert 1 AGN based on optical spectroscopy (|Duc et al.l 
[2003 ). This prior study of Abeh 1689 included mid-IR 
measurements from ISO CAM and the authors concluded 
that dusty star formation in this cluster is responsible for 
the vast majority of the observed mid-IR emission. Our re- 
sults support this conclusion as we find no other members 
with infrared AGN signatures. 

3.2.3. MS 1358+62 

Photometry an d spectroscop i c info rmation of MS 1358 
members is from IFisher et aLl () 19981 ). At z = 0.33, the 
3.3 and 6.2/tm PAH features shift into IRAC channels 2 
&; 4 respectively and colors along both axes are redder. 
This explains the dispersion seen in Figure [3] for late- type 
galaxies that are likely to be star forming. The spectro- 
scopic study of this cluster bv lFisher et al.l ()1998t ) revealed 
a number of emission line galaxies (ELGs). Not surpris- 
ingly, nearly all of the galaxies that depart from the "pas- 
sive cloud" are also ELGs. Of the galaxies in this cluster, 
we find one candidate AGN that is hosted by a late- type 
spiral. IMartini et all (|2009l ) find no X-ray AGN in this 
cluster with L(2-8 keV) > lO''^ ergs s~^. The one galaxy 
that we select as an AGN is detected as an ELG located 
roughly 860 kpc from the cluster center. 

3.2.4. CL 0024+17 

Spectrosco p y and photometry for CL 0024 is from 
IMoran et al.l (|2007a[) . star forming members with PAH 
emission are subject to the same effects as described for MS 
1358 and thus produce a similar scatter in mid-IR color- 
color space. The scatter here makes it difficult to discern 
the nature of the galaxies that are found near the bound- 
ary of the wedge. We do find two galaxies in the AGN- 
wedge (one appears to be sitting on the boundary). For the 
galaxy on the boundary, it is probable that star formation 



is producing its colors, though it is not ruled out as host- 
ing an IR-weak AGN. The other candidate AGN we find 
shows colors consistent with a po wer-law spectrum p lacing 
it well inside the IRAC wedge. I Zhang et all (|2005[ ) have 
conducted an X-ray observation of this cluster locating a 
handful of point sources at Lx > 10'*^ ergs s~^. However, 
none of the X-ray point sources overlap with any cluster 
galaxies from the spectroscopic catalog, i.e. these X-ray 
sources are in the field. 

CL 0024 is known to have numerous substructures 
as traced by thr ee techniques: weak-l ensing map 
(IKne ib et al.l 120031 ). X-ray shock fronts dZ hang et al.l 
2005) and a Dressler-Shectm an test ()Moran e t al. 2007at 
Dre ssier fc ShectmanI 119881] ). This implies that many 
galaxies are in groups that have been recently accreted (or 
are in the process) into the main cluster. Furthermore, the 
redshift distribution of galaxies shows a bimodality that 
suggests a rece nt merger along the line of sight with a large 
galaxy group (ICzoske et al.ll200ll 2002). The one candi- 
date AGN that is well inside the wedge has two close neigh- 
bors (confirmed spectroscopically) and is located '--^0.84 
Mpc from the cluster center (see Figure |4]). Redshifts of 
these galaxies are consistent with the main cluster. The 
candidate IR-AGN that we find near the edge of the IRAC- 
wedge is also within the main cluster; it is closer to the 
cluster core (~340 kpc) but is >80 kpc from the nearest 
neighboring galaxy. 

3.2.5. MS 0451-03 

Spectrosco py and photometry for MS 0451 are also from 
IMoran et aLl (l2]007a). At z = 0.54, the 6.2/im PAH feature 
shifts nearly outside the IRAC window. This reduces the 
scatter in mid-IR color due to star formation for members. 
X-ray data have shown that the distribu tion of the clus - 
ter ICM is predominantly elliptical (Donahue et al.l[2003l ). 
and the re dshift distribution o f cluster galaxies is broadly 
Gaussian (jMoran et al.ll2007al '). This indicates that MS 
0451 is predominantly virialized with no substantial in- 
falling galaxies. We detect no IR-AGN and no cluster 
galaxies are detected as X-ra y sources based on a Chandra 
survey (jMolnar et al.l [2002') , further indicating a lack of 
strong nuclear activity in MS 0451. 

3.2.6. MS 2053-04 

Spectroscopic and p hotometric informat ion of galaxies in 
this cluster come from lTran et al.l (|2005a[) . Detailed spec- 
troscopic and gravitational lensing stud ies of MS 2053-04 
([Verdugo et al.|[2007t ITran et al.|[2005aD show that it is a 
merger of two structures with 113 and 36 confirmed mem- 
bers respectively. Galaxies in the smaller structure (2053- 
B) have similar properties to field galaxies not associated 
with the cluster. This, coupled with the high fraction of 
star forming members (■~44%), suggests that MS 2053 has 
yet to completely virialize. 

Observations with Spitzer /MIPS also find that a fairly 
high fr action (~18%) of clus ter members are detected at 
24/im (jSaintonge et al.ll2008[ ). and the large population of 
star forming members clearly separates from the passive 
members in mid-IR color. Given the high level of activ- 
ity in MS 2053, it is surpising then that it only has one 
weak IR-AGN candidate. Using archival Chandra data, 
[Eastman et al.l (.2007.) find five X-ray sources at the red- 
shift of MS 2053 with Lx,h > 10^^ ergs/s. However, three 
are not classified as cluster members due to their distances 
from the cluster center (?■ > r2oo) and another is the BCG 
which is thought to be contaminated by X-ray emission 
from the ICM. This leaves only one cluster X-ray AGN 
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Fig. 4. — Thumbnails of our sample of infrared-selected AGN; images are from the Hubble Space Telescope. Numbers in parentheses 
correspond to the ID from Tablc[3] All images are ~150xl50 kpc and are oriented Nor th-up, East-left. Red circles correspond to the aperture 
of constant size (r 12.6 kpc) used to perform photometry on the IRAC images. See i|3.2.8l for an explanation of the labels in panel 8. 



which has a [3.6] — [4.5] color that is bluer than the IR- 
AGN selection region (Figure [3]). 



be identified as IR-AGN, it is not ruled out based on its 
[3.6] -[4.5] color (Figure E]). 



3.2.7. MS 1054-03 

Spectroscopica lly confirni e d mem bers and photometry 
are taken from ITran et aLl (|2007l ). One cluster X-ray 
source was not included in the photometric catalog, but 
optical and X-ray d ata for this galaxy are available from 
IMartel et al.l (|2007l ). Weak-lensing and X-ray analyses of 
MS 1054 sho w a clumpy natu re to the dark matter and 
ICM profiles (j.Tee et al.ll2005b[ ). The presence of such sub- 
structure indicates that the cluster experienced a merger 
and has yet to fully virialized. 

Earlier studies of MS 1054 reveal that it contai ns two 
members hosting X-ray AGN (j Johnson et al.l [20031 1 and 8 
radio sources that can be powered by AGN or star forma- 
tion (|Best et all [2001 ). We find only one infrared AGN 
at a distance >1 Mpc from the cluster center that is also 
detected as both an X-ray and radio source. The second 
X-ray source lies near the edge of the IRAC footprints and 
is not detected at 5.8/xm. Though this does not allow it to 



3.2.8. RX ,10152-13 



fror 



Spec troscopic mem b ership 



IBlakes 



of cluster galaxies is 
and photometry from 



Demarco etld] (|2005D , 

ee et al.l (|2006l) . This cluster shows signs of having 



gone through a large-scale merger event recently as 
indicated by its X-ray emi ssion, l uminosity distribution 
and weak-lensing profile (Je e et aL| [2005a'). 

We detect two IR-AGN in this cluster, both of which 
are classified as X-ray QSOs (Martel et al. 2007); no other 
cluster galaxies are detected as X-ray sources. One of 
the IR-K X-ray AGN is about 800 kpc from the cluster 
core while the second appears to be in a merging system. 
The latter detection is associated with five cluster galaxies 
within a projected radius of 30 kpc (Figure |4l neighbor- 
ing cluster galaxies are labeled). Due to their proximity 
and the IRAC PSF, these galaxies are blended into one 
mid-IR source where galaxy a is the closest (~0.7") to the 
centroid of the mid-IR emission. Redshifts in this quintet 
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TABLE 3 

IRAC-SELECTED CLUSTER, AGN 



ID 


Cluster 


RA (J2000) 


Dec (J2000) 


2^ 








HR° 


\og{LxT 


Morphology 


1 


Abell 1689 


13 11 35.5 


-01 20 12.8 


0.2000 


(1) 


0.30 


-21.66+0.04 






Sc 


2 


MS 1358 


13 59 24.0 


+62 31 08.0 


0.3236 


(2) 


0.86 


-21.64+0.03 






Sc 


3 


CL 0024 


00 26 40.0 


+17 09 41.8 


0.3955 


(3) 


0.33 


-21.36+0.05 






Sa+b 


4 


CL 0024 


00 26 33.7 


+17 12 19.8 


0.3964 


(3) 


0.84 


-21.64+0.05 






Sc+d 


5 


MS 2053 


20 56 21.0 


-04 37 22.8 


0.5763 


(4) 


0.19 


-21.48+0.10 






so/a 


6 


MS 1054 


10 57 02.7 


-03 39 43.6 


0.8319 


(5) 


1.02 


-24.12+0.07 


0.03+0.18 


43.23 


Irr 


7 


RX J0152 


01 52 43.8 


-13 59 01.3 


0.8201 


(5) 


0.78 


-23.61+0.08 


-0.62+0.05 


44.18 


Sb 


8 


RX J0152 


01 52 39.8 


-13 57 40.7 


0.8300 


(5) 


0.48 


-24.88+0.07 


-0.09+0.07 


44.52 


merger 


9 


RDCS 1252 


12 52 55.6 


-29 27 09.7 


1.2274 


(6) 


0.15 


-21.85+0.12 






Irr 


10 


RDCS 1252 


12 52 57.4 


-29 27 32.0 


1.2322 


(6) 


0.35 


-22.28+0.12 






Irr 


11 


RDCS 1252 


12 52 49.8 


-29 27 54.7 


1.2382 


(6) 


0.59 


-22.95+0.12 


0.20+0.31 


43.15 


Irr 


12 


RDCS 1252 


12 52 49.7 


-29 28 03.7 


1.2382 


(6) 


0.65 


-22.61+0.12 






SO 



Spectroscopic redshift: (1) Due et al. 2002 ; (2) Fisher et al. 1998 ; (3) Moran et al. 2007a ; (4) Tran et al. 2005a ; (5) Holden et al. 
2007 ; (6) Demarco et al. 2007. 

^ Projected distance from the cluster center in Mpc. 
Rest-frame 3.6/xm absolute magnitude. 

X-ray hardness ratio from lMartel et al.l | |2007I '| for AGN with X-ray detections. 
" Hard X-ray (2-10 keV) luminosity in ergs s~^ from Martel et al. 2007. 

f Morphology references: Michard in\ndrcon (2008, Coma), Due et al. (2002, Abell 1689l.'Fabricant et al.l ([200g . MS 1358), 'Moran et alj 
(20073, CL 0024), Moran et al. (2007a, MS 0451), Tr an et al.. (,,2005b . MS 20531. .Blakeslee et al. ( , 2006 . MS 1054) . iBlakeslee et al. (,2006, . 
RX J0152) and lDemarco et al.l (I2007l . RDCS J1252). 



are 0.867, 0.864, 0.834, 0.832 and 0.836 for galaxies a, b, c, 
d and e respectively. Based on these redsliifts, recessional 
velocities {w.r.t. to galaxy c) are 3743, 3408, 0, -231 and 
231 km/s for galaxies a, 6, c, d and e respectively. Since 
galaxies a k b have velocities within 400 km/s and galax- 
ies c, d e have velocities within 500 km/s (typical of 
galaxy groups) we suspect a h b are a bound system and 
c, d &: e are another bound system. Note, however, that a 
chance alignment of two groups such as this is '^0.1% likely 
to occur at random assuming a spherical cluster geometry 
with i? « 1 Mpc. Of course, due to the complex mor phol- 
ogy and substructure in RX J0152 (|Jee et al.ll2005al ) this 
probability may increase. 

3.2.9. RDCS J 1252-29 

At z = 1.24, RDCS 1252 is the most distant clus- 
ter in our sample yet has a virial mass similar to those 
of the lower redshi ft clusters as well as other properties 
(jRosati et al.l 120041) . Close inspection of the ICM in the 
cluster core reveals the presence of a shock front, sig- 
naling a recent merger with a cluster sub-clump. A de- 
tailed spectroscopic follow-up ([Demarco et al.l [ 20071) veri- 
fies a merger of two groups that have yet to viralize, and 
weak lensing shows that the centroid of the dark matter 
mass profile is offset from the optical/X-ray centroid by 
~8" (jLombardi et al.ll2005[) . Yet despite its young dynam- 
ical age, RDCS 1252 already has a population of luminous 
early-type galaxies that show little sign of ongoing star 
formation. 

At this redshift, emission from stellar populations be- 
gins to encroach into the IRAC window. This effect can 
be seen in Figure |3] as an upward shift in the the "passive 
cloud" relative to the lower- z clusters. This unfortunately 
brings the "passive cloud" closer to the AGN- wedge, possi- 
bly introducing contamination. Of the twenty nine RDCS 
1252 members shown in Figure[21 four are inside the AGN- 
wedge, of which one has been previously identified as an X- 
rayAGN (|Martel et al.|[2007KDe5iarco et al.ll2Q07[ ). Three 
of the four candidate IR-AGN are hosted by morpholog- 
ically irregular galaxies (Figure H]) that are likely to be 
gravitationally disrupted because of the large-scale clus- 



ter merger or by galax y-ga laxy mergers [HST im age from 
IBlakeslee et"air(|2003[ ) and IDemarco et al.l (|2003))- 

4. DISCUSSION 
4.1. Cluster IR-A GN Properties 

Recent studies of field galaxies find that infrared- 
selected AGN share similar properties ()Hickox et al.ll2009l : 
I Griffith fc Stern 2010): their host galaxies typically have 
late-type morphologies and tend to have blue optical col- 
ors. The properties of our 12 cluster IR-AGN agree with 
these earlier studies. We find the majority of our IR- 
selected AGN (10/12) are in late-type galaxies (Table [3]) 
that are blue (Fig. [5]) and so have recent/ongoing star for- 
mation. The remaining two cluster IR-AGN are hosted by 
SO galaxies that lie on/near the IR-AGN boundary (Fig.[2| 
and these members are consistent with being spiral galax- 
ies transi tioning to early-type systems in the clust er envi- 
ronment (jDressler et al.ill997l: iMoran et al.ll2007al) . 

It is worth noting that optical light from host galaxies 
may be contaminated by emission from an AGN, thus their 
colors could be biased bluewards. Our bluest AGN host, 
for example, is likely contaminated by the central engine 
since it is such a strong outlier. Figure S] shows ACS & 
WFPC-2 images of our IR-selected AGN, nearly all at rest- 
frame blue wavelengths. Upon careful visual inspection 
we find that the majority of these galaxies have extended 
morphologies and are not strongly dominated by a cen- 
tral point source. Furthermore. .Hickox et al.l (|2009D calcu- 
lated color contamination to be < 0.3 mag in °-^(m — r]3. 
When considering that the optical filters we use {B—V) are 
closer together in wavelength space and that AGN gener- 
ally contribute more flux at bluer wavelengths, we expect 
color contam ination in our sample to be less than what 
IHickox et al.l (HOOQ) find. Therefore, although point source 
contamination may be impacting some of our sample, we 
conclude that color contamination of our IR-AGN hosts is 
not significant enough to bias our results. 

* This color is computed by blueshifting the SDSS « & r filters by 
z = 0.1 
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Fig. 5. — Combined sky-plot for each cluster showing the projected 
locations of galaxies with respect to their cluster centers. IRAC- 
selected AGN are shown as circles (or stars for X-ray sources) and 
the dashrd fircle corresponds to a physical radius of 0.5 Mpc. Colors 
correspond to three bins is redshift, z < 0.5 (blue), 0.5 < z < 1.0 
(green) and z = 1.24 (red). A Kolmogorov-Smirnov test shows that 
the radial distribution of the IR-AGN is >99% likely to come from the 
same parent population as the normal cluster galaxies (gray points). 

We find that a third (4/12) o f our cluster IR- 
AGN are al so known X-ray sources IjMartel et all [2007t 
IJohnson et a l. 2003), a fraction that is nominally consis- 
tent with results from the AGN and Galaxy Evolution 
Survey ( AGES) where ~ 50 % of IR-AGN are also X-ray 
sources (jHickox et al.l [2009|) . However, it is worth men- 
tioning here that X-ray flux limits are not the same among 
these studies, and so agreement on this ratio is not neces- 
sarily implied. The X-ray sources are also the four most 
luminous 3.6/im IR-AGN (rest-frame; Figure |6]); however, 
note that the most luminous cluster IR-(- X-ray AGN is the 
blended source in RX J0152 (see §3.2.8). 

Figure [S] shows the projected distances of the 12 IR-AGN 
relative to confirmed members in all nine galaxy clusters. 
We find the radial distribution of the IR-AGN is drawn 
from the same parent population as the cluster galaxies 
with > 99% confidence using a Kolmogorov-Smirnov test. 
However, the four IR-I- X-ray AGN are all outside the clus- 
ter cores at Rproj ^ 0.5 Mpc. This result is consistent 
with 'E astman et~all (pOOTl) and lAtlee et al.l ([2QTTD who 
find that X-ray sources in galaxy clusters are not strongly 
centrally concentrated. Our observations may suggest that 
IR-I- X-ray AGN represent a different population than IR- 
only AGN; however, we are limited by the size of our sam- 
ple and so cannot further postulate on the uniqueness of 
these AGN based on their spatial distribution and 3.6/Ltm 
luminosities. 

4.2. Infrared- AGN Fractions 

To measure the fraction of IR-AGN in our cluster sample 
and test for evolution, we separate our sample into three 
redshift bins: low redshift (< 0.5), intermediate redshift 
(0.5 < z < 1.0), and a high redshift point at z = 1.24 
(RDCS 1252) containing 543, 377 and 29 IRAC detected 
members respectively. To ensure robustness, we consider 
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Fig. 6. — Rest-frame color-magnitude relation (optical color vs. 
3.6/im absolute magnitude); symbols are the same as in figure [5] 
Host galaxies of IR-AGN tend to have blue optical colors, thus these 
IR-AGN hosts have recent/ongoing star formation. Also, the four 
most luminous cluster IR-AGN are also X-ray sources. 



two different galaxy samples selected optically and in the 
mid-IR. We also take into account the varying spatial cov- 
erage of the IRAC mosaics and set the maximum field of 
view with the Coma cluster where the IRAC footprint only 
includes galaxies within Rproj ^ 760 kpc of the cluster cen- 
ter. In the higher redshift clusters, we therefore exclude 



members that are at Rproj > 760 kpc from their cluster 
center. 

Our first cluster galaxy sample is composed of optically- 
selected members brighter than V = —21.5 (this corre- 
sponds to where the ^-magnitude distribution turns over 
for RDCS 1252, our most distant cluster); this yields 118, 
141 and 22 galaxies in our three redshift bins. Note we 
do not to correct for passive evolution given that the host 
galaxies of the IR-AGN tend to have blue optical colors, 
i.e. are likely star forming systems. The IR-ACN fraction 
for this optically-selected sample is ~ 1% for both redshift 
bins at z < 1 and is only measurably non-zero at z = 1.24 
with /iR_AGN=18.2l^;'7% (Table H Fig. [71). AU errors in 
/iR,_AGN are asymni e tric l a Poisson uncertainties as de- 
termined bv lGehrelsl (|1986l ) for small number samples. 

Our second cluster galaxy sample is composed of mem- 
bers selected based on rest-frame 3.6/im luminosity. Be- 
cause of the evolution of the 3.6/im luminosity function 
is well-characteriz ed by passively evo lving galaxies that 
formed at z > 1.5 (jMuzzin et al.l[2008[ ). this luminosity se- 
lection is effectively a stellar mass cut. To isolate compara- 
ble samples of cluster members over our redshift range, we 
combine values of M*{z)[3.6iJ,m] (the characteris tic turn- 
ing point in the Schechter luminosity function; Schec hteil 
[1976)) from IMuzzin et all ([2008,) with the 80% limiting 
magnitude for our most distant cluster and thus select 
members brighter than rest-frame (Mg 5(2;) [3.6/im] -I- 0.5). 
We find that the cluster IR-AGN fraction is again uni- 
formly ~ 1% at 2 < 1 and only measurably non-zero at 
z = 1.24 with /iR-AGN=13.6;7^4% (RDCS 1252; Table H 
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Fig. 7. — Cluster IR-AGN fraction as a function of redshift for 
optically selected members brighter than Vab = —21.5 mag (green 
triangles) and mid-IR selected members brighter than (Mg g + 0.5) 
(blue circles). We consider three redshift bins: low redshift (z < 0.5), 
intermediate redshift (0.5 < z < 1.0), and a high redshift point at 
z = 1.24 (RDCS 1252) that contain 543, 377, and 29 IRAC-detected 
members respectively. The IR-AGN fraction is uniformly < 3% at 
2 < 1 and only measurably higher in RDCS 1252 a,t z = 1.24. Error 
bars represent Icr Poisson uncertainties derived using statistics from 
IGehrel: 
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Fig. ED. 

Thus far we have included all cluster galaxies regard- 
less of morphology in determining /ir_agn, but this 
may introduce a bias given that: 1) our cluster IR-AGN 
are predominantly hosted by late- type galaxies and 2) 
the morpholog i cal rn i x in clusters evolv e s with redshift 
(Dressier et al.' '19971: Fasano et all [2000t IPostman et"all 
,2005; Capak et al. 20Q3). In Fig.|8]we now exclude all mor- 
phologically classified E/SO members and measure a higher 
/iR-AGN at all redshifts (Tabled. However, /ir-agn re- 
mains < 5% at z < 1 in both of our selected galaxy sam- 
ples. Only in the most distant cluster (RDCS 1252) does 
/iR-AGN for late-typ^ members increase to ~ 70%. 

Although the number of cluster IR-AGN is small, we 
stress that our analysis is based on a sample of ~ 1500 
spectroscopically confirmed cluster galaxies at < z < 1.3, 
thus we place a strong upper limit on the IR-AGN fraction 
of < 3% for all members in massive clusters at z < 1. One 
caveat to consider is that while the IR color selection does 
identify - 90% of broad-line AGN it misses - 60% of 
narrow-line AGN (jStern et al.l 120051 ) and so we may be 
underestimating /ir-agn- However, strongly starbursting 
galaxies m ay also contarninat e our IR-AGN sample (up 
to ~ 50%: IDonlev et all 120081 ). and accounting for these 
actually decreases the IR-AGN fraction. Addressing these 
two competing effects is beyond the scope of our current 
analysis. 

The low IR-AGN fraction of - --^ 1% in ou r massive clus- 
ters at z < 1 is consistent with IMartini et al. (2009) who 
estimate an X-ray AGN fraction of 0.13 — 1.00% in clu sters 
at z 0.2 to z - 0.7 and with IHickox et al.l (|2009l ) who 

Here we mean all members except for E/SOs. 
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Fig. 8. — Same as Fig. [7] but only considering late-type galaxies 
(i.e. excluding E/SO galaxies). The IR-AGN fraction remains ~ 5% 
at z < 1 and is higher only at z = 1.24. 
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^ Members are selected using a luminosity limit in rest- frame 3.6/.tm 
or rest-frame Vab- 

^ Uncertainties in /ir-agn represent Icr Poisson errors determined 
from Gohrcls (1986). 

find the IR-AGN fraction in AGES is comparable to the X- 
ray AGN fraction at 0.25 < z < 0.8. However, we cannot 
say for certain that there is strong evolution in the cluster 
IR-AGN fraction with redshift given our small numbers. 
This is in contrast to the observed increase in the frac- 
tion of (d usty) sta r forming members in these same clus- 
ters (Saint onge et a l. 2008), thus the bulk of their 24/xm 
flux is due to star formation and not AGN. Our single 
galaxy cluster at z > 1 does suggest that IR-AGN have a 
more prominent role at this epoch, but we recognize that 
1) RDCS 1252 may be unusually active and 2) the IRAC 
color selection starts to suffer contamination from passive 
galaxies at these redshifts. Given their rarity, a larger sur- 
vey of IR-AGN in massive galaxy clusters, particularly at 
z > 1, is needed to robustly identify any evolution in IR- 
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AGN with redshift. 

Another interesting comparison we can study is the vari- 
ation in /iR-AGN with environment. Using magnitude cuts 
similar to our V < —21.5 hmit, the /ir-agn in the Bootes 
field sample from the AGES (0.25 < z < 0.8) is -2% (R. 
Hickox, private communication). This is well within our 
upper Icr uncertainty (Table H]) at similar redshifts, thus 
there is no clear variation in /ir_agn based on local den- 
sity. The possibility of AGN playing a more influential role 
at 2 > 1 still remains. 



5. CONCLUSIONS 

We present the first census of mid-infrared selected ac- 
tive galactic nuclei (IR-AGN) in massive galaxy clusters 
(A^vir > 5 X IO^Mq) at < z < 1.3 by combin- 
ing archival Spitzer /IRAC imaging with extensive optical 
spectroscopic catalogs (public and private) and deep op- 
tical photometry of ~ 1500 confirmed members in nine 
clusters. Our clusters are selected to be the most mas- 
sive well-studied systems currently known. Using the 
four IRAC channels (3.6, 4.5, 5.8 and S .O^m) and estab- 
lished mid-IR col or selection techniques (jStern et al.1120051 : 
ILacv et aLl I2004D . we identify 949 members that are de- 
tected (> 3(t) in at least three of the four IRAC channels 
and isolate 12 that host dust-enshrouded AGN. Similar 
to IR-selected AGN in rece nt field studies (|Hickox et al.l 
[2009HGrifRth fc S tern 20100. the host cluster galaxies tend 
to be late-type members with blue optical colors that in- 
dicate recent/ongoing star formation. The IR-AGN have 
the same radial distribution as the cluster members, but 
the four most IR-luminous AGN lie outside of their clus- 
ter cores {Rproj > 0.5 Mpc) and are also known X-ray 
sources. This suggests that very bright IR-I- X-ray AGN 
are not centrally concentrated in their clusters, consistent 
with t he results for bright X-ray sources by Martin Fet al.l 
(|2009f) . Our results suggest that IR-|-X-ray AGN may not 
be the same population as the IR-only AGN, but we are 
too limited by our sample to make any assertion. 

To measure the fraction of IR-AGN and test for evolu- 
tion, we compare two complete samples of cluster galax- 
ies: 1) an optically-selected sample with members brighter 
than Vab = —21.5 (rest-frame) and 2) a mid-IR selected 
sample with members brighter than (M*(z)[3.6/im] -I- 0.5) 
(jMuzzin et a l. 2008) that is essentially a stellar mass cut. 
For the eight galaxy clusters at 2 < 1, we place a strong 
upper limit of < 3% on the fraction of IR-AGN for both 
cluster samples. Because IR-AGN tend to be hosted by 
late-typ e galaxies and the morphological mix in clusters 
evolves (jPressler et al.l[T997l : IPostman et al.l[2005[ ) , we also 
consider only late-type members and find that the frac- 
tion with IR-AGN is < 5% for both samples. These low 
IR-AGN fractions are surprising given that the fraction 
of (dusty) star formation in these same c lusters increases 
by about a fact or of four at < z < 1 (jSaintonge et al.l 
[200a iBaTeFan [2010). However, an IR-AGN fraction of 
~ 1% is consistent with the low fr action of X-ray AGN in 
galaxy clusters (< 1% M artini et a l. 2 009.) and the rela tive 
populations of X-ray vs. IR AGN t Hickox et al.|[2009[ ). 

In contrast, our single galaxy cluster at z = 1.24 (RDCS 



1252) has a measurably higher IR-AGN fraction of ^ 15% 
(all galaxy types) and ~ 70% (late-types only). However, 
RDCS 1252 may simply be an unusually active cluster. 
Also, the IR color selection starts to suffer stronger con- 
tamination from non-AGN members at z > 1.2. 

We also compare our /ir-agn measurements in dense 
clusters at z < 1 to that of the Bootes field from the AGES, 
which probes sparser galactic environments over a similar 
range in redshift (0.25 < z < 0.8). Using optical magni- 
tude cuts similar to this study, /ir-agn is measured to 
be '~2% in the field (R. Hickox, private communication). 
Consequently, we do not see a statistically significant vari- 
ation in /iR_AGN here that would be caused by local galaxy 
density. However, the question still remains as to whether 
or not IR-AGN have a more profound impact at z > 1. 

We note that while the IR color selection successfully 
identifies ~ 90% of b road-line AGN, it does miss ~ 60% 
of narrow-line AGN (jStern et al.l 120051 ) and so /ir-agn 
is undoubtedly incomplete, i.e. underestimated. On the 
other hand, contamination may arise from star forming 
galaxies falsely identified as IR-AGN; thus we may also 
be overestimating /ir_agn- Such co i itamination can be 
as high as 20-50% (IDonlev et al.l[2008l : iHlckox et al.l[2009[ ) 
and is more influential at low luminosities. For simplicity, 
we ignore these two competing effects because they are 
beyond the scope of our analysis and do not change our 
general conclusions. 

Taken as a whole, our results show that IR-AGN and 
star formation are not strongly correlated at z < 1 be- 
cause the IR-AGN fraction is uniformly very low (~ 1%) 
at z < 1 whereas sev eral of these clusters h ave star forming 
fractions of > 20% (|Saintonge et al.l[2008[) . We do find a 
hint of evolution in the IR-AGN fraction at z ~ 1.2, but 
only with a more extensive mid-IR survey of galaxy clus- 
ters, particularly at z > 1, can we confirm this intriguing 
result. A study with equally deep spectroscopic coverage 
across this wide redshift range is also needed to robustly 
measure /ir_agn and test for evolution as a function of 
environment. While Spitzer /IRAC is no longer available, 
such studies should be possible with the upcoming public 
survey by the Wide-field Infrared Survey Explorer (WISE) . 
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